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Abstract 

A square matrix is nonderogatory if its Jordan blocks have distinct 
eigenvalues. We give canonical forms for 

• nonderogatory complex matrices up to unitary similarity, and 

• pairs of complex matrices up to similarity, in which one matrix 
has distinct eigenvalues. 

The types of these canonical forms are given by undirected and, re- 
spectively, directed graphs with no undirected cycles. 
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1 Introduction 

A square matrix is nonderogatory if its Jordan blocks have distinct eigenval- 
ues; that is, if its characteristic and minimal polynomials coincide. 
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We give canonical forms for 

• nonderogatory matrices up to unitary similarity, and 

• pairs of matrices up to similarity, in which one matrix has distinct 
eigenvalues. 

All matrices that we consider are complex matrices. 

Our canonical matrices are special cases of the canonical matrices that 
were algorithmically constructed by Littlewood and Belitskii: 

• Littlewood [10] developed an algorithm that reduces each square matrix 
M by unitary similarity transformations 

is a unitary matrix, 

to a matrix Mean in such a way that M and are unitarily similar if 
and only if they are reduced to the same matrix Mean = A^can- Thus, 
the matrices that are not changed by Littlewood's algorithm are canon- 
ical with respect to unitary similarity. Other versions of Littlewood's 
algorithm were given in |1] and [121 E] • 

• Belitskii [H |2] developed an algorithm that reduces each pair of n x n 
matrices (M, A^) by similarity transformations 

(M, N) ^ {S-^MS, S-^NS), S is nonsingular, (1) 

to a matrix pair (M, N)can in such a way that (M, A^) and (M', A^') 
are similar if and only if they are reduced to the same matrix pair 
(M, A^)can = (M', A^')can- Thus, the matrix pairs that are not changed 
by Belitskii's algorithm are canonical with respect to similarity. Be- 
litskii's algorithm was extended in [TS] to the problem of classifying 
arbitrary systems of linear mappings and the problem of classifying 
representations of finite dimensional algebras. 

Lists of Littlewood's canonical 5x5 matrices and Belitskii's canonical 
pairs for 4 X 4 matrices are in [H] and Without restrictions on the size 
of matrices, we cannot expect to have explicit descriptions of Littlewood's 
canonical matrices and Belitskii's canonical matrix pairs since 

• The problem of classifying matrices up to unitary similarity contains 
the problem of classifying arbitrary systems of linear mappings on uni- 
tary spaces [Ulllj; and 
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• The problem of classifying matrix pairs up to similarity contains the 
problem of classifying arbitrary systems of linear mappings on vector 
spaces [HI |3] . 

When it is applied to nonderogatory matrices, Littlewood's algorithm can 
be greatly simplified. Mitchell [TT] presented an algorithm intended to reduce 
nonderogatory matrices to canonical form, but his algorithm is incorrec t0. In 
Sections |2H5] we give a version of Littlewood's algorithm for nonderogatory 
matrices and describe a set of canonical nonderogatory matrices for unitary 
similarity. Each type of canonical nonderogatory matrices with t distinct 
eigenvalues is given by an undirected graph with t vertices and no cycles. 

When it is applied to pairs of n x n matrices in which one matrix has 
distinct eigenvalues, Belitskii's algorithm can also be greatly simplified. In 
Section [6] we describe a set of canonical forms for pairs of matrices under 
similarity. It is analogous to the set of canonical nonderogatory matrices in 
Section HI but it involves directed graphs instead of undirected graphs. This 
description was used in [5] to classify pairs of 4 x 4 matrices up to similarity. 

2 Schur's triangular form for nonderogatory 
matrices 

Schur's unitary triangularization theorem [71 Theorem 2.3.1] ensures that 
each square matrix M is unitarily similar to an upper triangular matrix 

Ai ai2 . . . ain 

'^n— 1,71 



^The following reasoning on page 71 of P~l is incorrect: "Let us agree to go from left to 
right down the successive diagonals below the main diagonal and pick out each non-zero 
element as we come to it until we obtain either a total of n — 1 non-zero elements or all non- 
zero elements off the main diagonal, where n is the order of the matrix. These chosen non- 
zero elements can then be made positive by transforming by a diagonal unitary matrix." 
Unfortunately, it is impossible to make positive in this way "each non-zero element as we 
come to it" . In Section [3] we choose a set of nonzero elements that can be made positive. 



Ai < As < • • • < A„, (2) 
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whose diagonal entries are complex numbers in any prescribed order; for 
definiteness, we use the lexicographic order: 



a + bi ^ c + di if either a < c, or a = c and b ^ d. 



(3) 



A unitary matrix U that transforms M to an upper triangular matrix 
A = U^^MU of the form (|2]) can be constructed as follows: first find a 
nonsingular matrix S such that J = S^^MS is the Jordan form of M that 
has diagonal entries in the prescribed order, then apply the Gram-Schmidt 
orthogonalization to the columns of S and obtain a unitary matrix U = ST, 
in which T is upper triangular. Alternatively, a unitary U with the desired 
property can be constructed directly, without first obtaining the Jordan form 
[TJ Theorem 2.3.1]. 

The unitary similarity class of M can contain more than one upper trian- 
gular matrix A of the form For example, the argument of any nonzero 
entry in the first superdiagonal may be chosen arbitrarily. The following 
diagonal unitary similarity permits us to standardize the choice of these ar- 
guments by replacing every nonzero entry tti^i+i in the first superdiagonal by 
the nonnegative real number rj := |aj^j+i|: 



A ^ UAU 



-1 



U := diag(l, Ui, UiU2, U1U2U3, 



in which Ui 



^1,1+ 



i/vi if tti^i+i ^ and Ui := 1 if ai^i+i = 0. This unitary 



similarity is used in the following example. 

Example 2.1. Every square matrix M that is unitarily similar to a matrix of 
the form ([2]), in which all entries of the first superdiagonal of A are nonzero, 
is unitarily similar to a matrix of the form 



B 



Ai 612 
A2 





'In 



bn—l,n 
An, 



all bi^i+i > 0, 
Ai ^ A2 ^ • • • ^ A„. 



(4) 



Such a matrix can be used as a canonical form for M under unitary simi- 
larity since if two matrices of the form (|4]) are unitarily similar, then they 
are identical. This canonical form is a special case of a canonical form for 
nonderogatory matrices that we construct in Section |H 
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The number of Jordan blocks with eigenvalue A in the Jordan form of an 
n X n matrix A is equal to n — rank(A — A/„). Thus, a matrix of the form 
([2]) is nonderogatory if and only if Aj = Aj+i implies that aj^j+i 7^ 0. We 
formalize this observation in the following lemma. 

Lemma 2.1. A matrix is nonderogatory if and only if it is unitarily similar 
to a block matrix of the form 



A 



Ai A12 
A2 





A 



It 



At-i,t 
A* 



(5) 



in which each diagonal block Aj is nii x rrii and has the form 
\j * ... * 



A,- 



A,: 







* 

A. 



all entries of the first superdiagonal 
of Ai are positive real numbers, 



(6) 



and the diagonal entries are lexicographically ordered: Xi < X2 < ■ ■ ■ < Xt- 



3 An algorithm for nonderogatory matrices 



Let M be a nonderogatory matrix. We first reduce it by unitary similarity 
transformations to a matrix A of the form described in Lemma 12.11 Then we 
reduce A by transformations A^ A := U^^AU {U is unitary) that preserve 
this form. 

We prove in Lemma 15.11 that A' has the form described in Lemma 12.11 if 
and only if 



Thus, we reduce A to canonical form by transformations 

'Ai u^^U2Ai2 . . . u^'^utAit 

A2 ■■• ; 

■ • . u~\utAt-i^t 



1. 



A 







A, 



(7) 
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Notice that the blocks Aij are multiphed by complex numbers of modulus 1. 

We construct a set of canonical nonderogatory matrices that includes the 
canonical matrices from Example 12.11 For this purpose, if A12 0, then we 
reduce it to the following form. 

Lemma 3.1. Let C = [cij] be a nonzero p x q matrix. Let c be the first 
nonzero entry in the sequence formed by the diagonals of C starting from the 
lower left: 

Cplj Cp— 1^1, Cp2i Cp— 2,1) Cp— 1^2) Cp3, . . . , C\q. 

We can replace c by the positive real number r = |c| by multiplying C by a 
complex number of modulus 1 . The resulting matrix is canonical with respect 
to multiplication by complex numbers of modulus 1. 

For example, if the first nonzero diagonal (starting from the lower left) 
of C is below the main diagonal, then its canonical matrix from Lemma [3.11 
has the form 



* 

■•• 
■•• * 
■•. 

Or* 

* * 

_ * * 

We sequentially reduce the blocks A^j of the matrix ([5]) to canonical form 
in the following order (i.e., arranging them along the block superdiagonals of 
A): 

A12, A23, At-i^t'-i Ai3, A23, At-2,t', Alt. (8) 

We begin with the block A12. If A12 = 0, then it is not changed by 
transformations of the form (JTj), and so it is already canonical. If A12 0, 
then we reduce it as in Lemma 13. Ij to preserve the block A12 obtained, we 
must impose the condition Ui = U2 on the transformations ([7]). 

Then we reduce ^23 in the same way and so on, until all blocks in the 
first superdiagonal have been reduced. We obtain a matrix A in which all 
nonzero blocks in the first superdiagonal have the form described in Lemma 



r G M, r > 0, 
*'s are complex numbers. 
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13.11 This matrix is uniquely determined by the unitary similarity class of A, 
up to transformations of the form ([7]) that satisfy the conditions Ui = Ui+i if 
^ 0; we say that such transformations are admissible. It is convenient 
to describe these conditions by a graph G*-^^ with vertices 1, . . . ,t and with 
edges i — + 1) that correspond to all A^j+i ^ 0. 

Next we reduce the blocks of the second superdiagonal to canonical form. 
If = or if lii = ^2 = Ms (i.e., if G'-^^ contains the path 1 — 2 — 3), then 
^13 is not changed by admissible transformations of the form ([7j); it is already 
canonical. If Ais ^ and G^^^ does not contain the path 1 — 2 — 3, then 
we reduce as in Lemma 13.11 and add the edge 1 — 3 to the graph. Then 
we reduce A24 and so on until we have reduced all blocks in the sequence ([H]). 

This algorithm can be formalized as follows. For each graph G with 
vertices 1, . . . , t, we say that ([7j) is a G -transformation if Ui = uj for all edges 
i — j in G. 

Algorithm 3.1. Let M be a nonderogatory matrix, let A be its upper 
triangular form ([5]) for unitary similarity described in Lemma [2.11 and 
let Go be the graph with vertices 1, . . . ,t and without edges. 

The first step: We construct a pair {Ai, Gi) as follows. Let A^^^^ be the first 
nonzero block of A in the sequence ([8]). Reduce Ap^qi as in Lemma [3TT] 
by transformations of the form ([7]) and denote the resulting matrix by 
Ai. Add the edge pi — qi to Gq and denote the resulting graph by Gi. 

The ath step (a ^ 2): Using the pair {Aa-i,Ga-i) constructed at the (a — 
l)st step, we construct {Aa,Ga)- Let Ap^g^ be the first block of 
Aa-i that is to the right of in (jS]) and is changed by G^-i- 

transformations (this means that Ap^g^ 7^ and Gq,_i does not contain 
a path from pa to Qa)- We reduce Ap^g^ as in Lemma [3?T] and denote 
the resulting matrix by Aa- Add the edge pa — Qa to Gq_i and denote 
the resulting graph by Gq,. 

The result: The process stops at a pair (^4^, Gr) such that all blocks of Ar to 
the right of Ap^g^ in ([8]) are not changed by G^-transformations. The 
number r of steps is less than t since the graph Gr has t vertices, r 
edges, and no cycles. Write Mean := A^. and G := G^. 

In the proof of Theorem 14. II we show that the pair (Mean, G) is uniquely 
determined by the unitary similarity class of M; that is. Mean is a canonical 
form for M with respect to unitary similarity. 
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4 Canonical nonderogatory matrices and the 
classification theorem 



Algorithm 13. II constructs a pair (Mean, G) for each nonderogatory matrix M. 
The structure of Mean is determined by the graph G as follows: 

• The blocks 

of Mean have the form described in Lemma I3.1j they correspond to the 
edges of G. 

• Let Aij [i < j) be a block of Mean that is not a member of the list in 
OH]). Let Ap^g^ be the nearest block in the list that is to the left of 
Aij in dH]). If there is no such block (i.e., if Aij is to the left of Ap-^q-^), 
we put a := 0. Then 

(i) Aij = if Ga does not contain a path from i to j, and 

(ii) Aij is arbitrary if Ga contains a path from i to j. 

The graph G^ in (i) and (ii) can be obtained from G by removing the edges 
u — V that correspond to those Auv in the list that are reduced after 
Ap^q^ if a 7^ 0, and by removing all the edges of G if a = 0. Thus, Auv is to 
the right of Aj in dH]); i-e., either v — u>j — i, or v — u = j — i and u > i. 
Hence, Mean is a G-canonical matrix in the sense of the following definition. 

Definition 4.1. Let G be an undirected graph with vertices 1,2, ... ,t and 
no cycles. By a G-canonical matrix, we mean a block matrix of the form ([5]) 
in which every diagonal block has the form ([6]), Ai < A2 < • • • < At, and each 
block Aij {i < j) satisfies the following conditions: 

• Aij has the form described in Lemma [3. II if G contains the edge i — j; 

• Aij = if either G contains no path from i to j, or the path from i 
to j (which is unique since G without cycles) contains an edge u — v 
{u < v) such that 

— either v — u > j — i, 

— or v~u = j — i and u > i; 
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• Aij is arbitrary, otherwise. 
Example 4.1. (a) Each matrix of the form (jl]) is G-canonical with 

G : 1 2 1 

(b) Each G-canonical matrix with 

5 

G : 

2 3 4 1 

has the form 





in which 

— each block Aj has the form ([6]) and Ai < A2 < A3 < A4 < A5, 

— each block Cj has the form described in Lemma [3.11 

— the stars denote arbitrary blocks. 

A G-canonical matrix is a canonical nonderogatory block if G is a tree. 
It follows from the uniqueness in (b) of the next theorem that canonical 
nonderogatory blocks are indecomposable under unitary similarity, i.e., they 
are not unitarily similar to a direct sum of square matrices of smaller sizes. 
Their role is analogous to the role of Jordan blocks in the Jordan canonical 
form. 

Theorem 4.1. (a) For each nonderogatory matrix M, there is a unique 
undirected graph G and a unique G-canonical matrix Mean such that M is 
unitarily similar to Mean- Thus, M is unitarily similar to N if and only if 

■Mean A'can • 

(b) Each nonderogatory matrix M is unitarily similar to a direct sum of 
canonical nonderogatory blocks. This direct sum is uniquely determined by 
M , up to permutation of summands. 









C4 




A2 




* 


* 




A3 


C2 


C3 






A4 











As 
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Remark 4.1. The direct sum in Theorem 14.1( b) is permutationally similar 
to Mean and can be obtained from it as follows: The graph G is a disjoint 
union of trees; denote them by Gi, . . . ,Gs. Let un < Ui2 < ■ ■ ■ < uu- be 
the vertices of Gi. Let Ai be the ti x ti sub matrix of Mean formed by rows 
Uii, . . . , Uiti and columns un, . . . , uu^. Definition 14.11 ensures that the uu.Ujk 
block of Mean IS zero if z 7^ j. Therefore, Mean is permutationally similar to 

Ai®A2®---®A,, 

which is the desired direct sum. Each Ai is a ^--canonical matrix, in which 
G[ is the tree obtained from Gi by relabeling the vertices un, . . . ,Uit- with 
1, . . . , tj. 



5 Proof of Theorem 4.1 



Our proof is based on the following lemma about unitary similarity of matri- 
ces of the form described in Lemma 12. 1[ This lemma was proved in greater 
generality in [10] and in p4]; we offer a proof for the reader's convenience. 



Lemma 5.1. Let 

"Ai Ai2 

A = 





A 



It 



A-2 



At-u 
At . 



B 



a; 

a; 



Bit' 







Bf-i^t' 

K' . 



he nonderogatory matrices of the form described in Lemma lKTi Assume that 
they are unitarily similar: U^^AU = B with unitary U . Then t = t' , 



Ai = a;, 



A, = A' 



and U has the form 



U = Ml/mi © • • • © Utl„ 



(10) 



(11) 



in which ui, . . . ,ut are complex numbers of modulus 1 and the size of Aj is 
^ T^i for each i. 

Proof. The matrices A and B have the same main diagonal since they are 
similar and the entries along their main diagonals are lexicographically or- 
dered. This means that t = t' and for each i the diagonal blocks Aj and A^ 
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are rrii x rui matrices of the form (E]) with the same Aj 
into three steps. 

Step 1: Prove that U has the form 

U =Ui®U2®---®Ut 



The proof is divided 



(12) 



in which every block Ui is rUi x nii. If t = 1 there is nothing to prove, so 
assume that t ^ 2. Partition U into blocks Uij of size mj x rrij. Our strategy 
is to exploit the equality of corresponding blocks of both sides of the identity 
AU = UB. 

The t,l block of AU is AtUn, and the t,l block oi UB is UtiA[. Since 
At 7^ Ai, Uti = is the only solution to AtUn = UtiA[. 

If t > 2, then the t,2 block of AU is AtUt2, and the t,2 block of UB 
is Ut2A'2 (since Ua = 0); we have AtUt2 = Ut2A'2- Since \t A2, we have 
Ut2 = 0. Proceeding in this way across the last block row of AU = UB, we 
find that Uti, Ut2 ■ • • , Ut,t-i are all zero. 

Now equate the blocks of AU = UB in positions {t — l),k for k = 
1, 2, . . . , t — 2 and conclude in the same way that t/t-1,21 ■ ■ ■ , ^t-i,t-2 

are all zero. Working our way up the block rows of AU = UB, left to right, 



we conclude that Uj = for all i > j. Since U 



-1 



U*, it follows that 



for all j > i and hence 



U =Uu®U22®---®Uu. 



This proves ( IT^ with U := Uu. 

Step 2: Prove that U is diagonal. Since AU = UB, we have t identities 
AiU = UiA[, i = 1, . . . ,t, and all the entries in the first superdiagonal of 
each Aj and are positive real numbers. Thus, it suffices to consider the 
case t = 1. In this case 



A 



A ai2 
A 





Q-n— 1,71 



A 



B 



A 612 
A 





bin 



-'n— l,n 

A 



tti^i+i and hi^i+i are positive real numbers for all z = 1, . . . , n — 1, and AU 
UB. As in Step 1, we equate corresponding entries of the identity 



{A-\In)U =U{B~\In 



(13) 
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In position n,l we have = 0. In position n,2 we have = 'Uni^i2; since 
612 7^ it follows that Uni = 0. Proceeding across the last row of (IT^ . we 
obtain 

Unl = Un2 = ■ ■ ■ = Mn,n-1 = 0. 

Working our way up the rows of f|T3|) in this fashion, left to right, we find 
that Uij = for all i > j. Thus, U is upper triangular. Since U is unitary, it 
is diagonal: U = diag(ui, . . . , 

Step 3: Prove that U = diag(Mi, . . . , Un) has the form (fTT]) . We continue 
to assume that t = 1. Equating the entries of AU = UB m position i,i + 1, 
we have aj^j+iitj+i = Uibi^i+i, so aj^j+i/6j^j+i = Ui/ui+i, which is positive real 
and has modulus one. We conclude that Ui/ui+i = 1 for each z = 1, . . . , n— 1, 
and hence «! = ••• = «„. This proves fllip . which implies (ITU]) . □ 

Proof of Theorem \4.1\ (a) Let M be a nonderogatory matrix. Algorithm 13. II 
constructs the graph G and the matrix Mean, which is unitarily similar to 
M. As shown at the beginning of Section HJ Mean is a G-canonical matrix. 

Let M and be nonderogatory matrices that are unitarily similar. Our 
goal is to prove that Algorithm 13 . 1 1 reduces them to the same matrix Mean = 
A'ean- Followiug the algorithm, we first reduce M and N to matrices A and 
B of the form described in Lemma 12.11 They are unitarily similar; that is, 
U^^AU = B for a unitary matrix U . Lemma [5. II ensures that t = t', Aj = A^ 
for all i, and U has the form f lTT]) . This means that B is obtained from A by 
a transformation of the form ([71): 

Bij = u^^UjAij, \ui\ = 1, i,j = l,...,t. (14) 

We arrange the superdiagonal blocks Aij in A and Bij in B along the 
block superdiagonals, as in ([S]). By f[T^ . the first nonzero superdiagonal 
block of A and the first nonzero superdiagonal block of B occur in the same 
position pi,qi. In Step 1 of Algorithm 13.1] we reduce them to the same form 
described in Lemma l3.ll and obtain the matrices Ai and Bi, in which the 
Pi,qi blocks are equal. 

In Step a, we reduce the first superdiagonal block of Aa-i that is changed 
by Go-i-transformations, and the first superdiagonal block of Ba-i that is 
changed by Ga„i-transformations. They occur in the same position Pa^Qa 
and are reduced to the same form described in Lemma 13.11 We obtain the 
matrices A^ and B^, in which the blocks in position Pa^Qa coincide; the 
superdiagonal blocks that precede them coincide as well. The matrix B^ can 
be obtained from A^ by a Go-transformation, which preserves these blocks. 
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The process stops at a matrix such that none of its blocks are changed 
by Gr-transformations. Then Ar = and so M^an = -^can- 

(b) This statement follows from Remark 14.11 □ 



6 Canonical matrix pairs for similarity 



Let (M , A^) be a pair of nxn matrices, and let M have n distinct eigenvalues. 
In this section, we give a canonical form for (M, A^) with respect to the 
similarity transformations ([1]). 

The pair (M, A^) is similar to some pair 



( 



\ 



Ai 











A. 



&11 



bir 



\ 



J 



Ai < ■ • ■ < A, 



(15) 



in which < is the strict lexicographic order on C; see 
Let 



{A', B') 



a; 











A' 



b' 



In 



b'nn 



a'i<---<a;, 



be another pair of this form, and let it be similar to {A, B); that is, 
{S-^AS,S-^BS) = iA',B') for some nonsingular S. Then A = A', 
AS = SA, and so 5" = diag(si, . . . , s„) in which e C. Thus, 

the pair f lT5]) is uniquely determined by {M,N), up to transformations 



B^ B' 



bii ^S2bi2 ... ^Snbin 

S2^Sib21 622 • • • S2^Snb2n 

Sn^Sibnl S^^S2bn2 ■■■ bnn 



(16) 



in which Si, . . . ,Sn are arbitrary nonzero complex numbers. 

Example 6.1. Suppose that a pair (M, N) of n x n matrices is similar to a 
pair of the form ( ITSl) in which 612, 613, ... , 6in are all nonzero. Taking si = 1, 
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S2 = bi2, ■ ■ ■ , Sn = b^n dlS]), we reduce (M, A^) to the form 



V 



Ai 







A, 



Xr. 



* 1 ... 1 

* * ... * 



Ai < • • • < A, 



(17) 



in which the stars denote arbitrary complex numbers. We can use f[T7j) as a 
canonical form for (M, N) for similarity since if B and B' in flTB]) have 1 in 
positions l,k, k = 2,3, ... ,n, then Si = ■ ■ ■ = and so i? = B'. Thus, if 
pairs of the form (ITTI) are similar, then they are equal. 



In the general case, we reduce B by transformations of the form flTBl) 
using the following algorithm. We arrange the entries of B along the rows 
starting from the first; that is, bij precedes bpq if {i,j) < {p,q) with respect 
to the lexicographic order. For each directed graph G with vertices 1, . . . , n, 
we say that (fT6!) is a G -transformation if Si = sj for all directed edges i ^ j 
in G. 

Algorithm 6.1. Let B = [bij] be an t?, x -n, matrix. Denote by Go the graph 
with vertices 1, . . . ,n and without edges. 

The first step. The entry bu is not changed by transformations of the form 
(!T6|) : we mark it as reduced and write {Bi, Gi) := {B, Go). 



The second step. If 612 = then it is not changed by Gi-transformations, we 
mark 612 as reduced and write {B2, G2) := (-Bi, Gi). If 612 7^ then we 
make 612 = 1 by Gi-transformations, add the directed edge 1 ^ 2 to 
Gi, and denote by B2 and G2 the resulting matrix and directed graph. 

The kth step. Let bpq be the kth entry; that is, {p — l)n + q = k. li p = q, 
or bpq = 0, or G^-i has an undirected path from p to q, then bpq 
is not changed by Gfc_i-transformations; we mark bpq as reduced and 
write {Bk,Gk) ■= {Bk~i,Gk~i). Otherwise, we make bpq = 1 by Gu-i- 
transformations, add the directed edge p — * g to Gfc_i, and denote by 
Bk and Gk the resulting matrix and directed graph. 

The result. After steps, we obtain a matrix i?„2, in which all entries have 
been marked as reduced. Write (i?can,G) := (S„2,G„2). 
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Let us show that -Bean ^-5 a canonical form for B with respect to transfor- 
mations of the form f ll6p : that is, if B and C are n x n matrices such that B 
can be reduced to C by transformations of the form f|T6l) then Scan = C'can- 
Indeed, after k steps of Algorithm 16.11 apphed to B and C, we obtain the 
matrices Bk and Ck and the same directed graph Gk- One can prove by 
induction on k that Bk reduces to Ck by G^-transformations, and so the first 
k entries of Bk and Ck coincide. Taking k = n"^, we obtain B^^n = C'can- 

Let (M, A^) be a pair of n x n matrices, and let M have n distinct eigen- 
values. Then {M,N) is similar to a pair {A,B) of the form (fT5|) . which 
is uniquely determined by {M,N), up to transformations of the form (fT6|l . 
Taking S = -Beam we obtain the pair (M, A)can := (^i -Bean), which is similar 
to {M,N) and is uniquely determined by {M,N). Thus, 



(M, A)can a canonical form for (M, A^) /or similarity. 



(18) 



In the fcth step of Algorithm |6lTl we reduce the kth entry bpg and construct 

by 
that 



the directed graph Gk- The graph Gk can be also obtained from G 
removing the directed edges i j that correspond to those entries bij 



Gn2 



were reduced to 1 after bpg] this means that 



,j)> {p,q). Thus, 



the pair {M,N)can is G -canonical 
in the sense of the following definition. 



(19) 



Definition 6.1. Let G be a directed graph with vertices 1,2, ... ,n and no 
undirected cycles. By a G-canonical matrix pair we mean a matrix pair of 
the form 



Ai 











A. 



foil 



Ki 



bir 

brr 



Xl ^ • • • ^ Xj., 



in which every entry bpg satisfies the following conditions: 
fi) b„n = 1 if G has the directed edge p ^ q; 



n 



^pq 



if either G has no undirected path from p to q, or the undirected 
path from ptoq contains a directed edge i ^ j such that {i,j) > {p, q) 
with respect to the lexicographic order; 

iii) bpg is arbitrary, otherwise. 



15 



Example 6.2. (a) Each pair of the form fll7l) is G-canonical with 




(b) Each G-canonical matrix pair with 



G : 



has the form 



v 



Ai 

















As 

















A3 

















A4 

















As 



* 1 
1**0 

* * * 
1* 



in which Ai < ■ ■ ■ < A5 and the stars denote arbitrary complex numbers. 

A G-canonical matrix pair is an indecomposable canonical matrix pair if 
G is a tree. It is not similar to a direct sum of pairs of square matrices of 
smaller sizes. This is a consequence of the uniqueness assertion in (b) of the 
following theorem. 

Theorem 6.1. (a) For each pair {M,N) of n x n matrices in which M has 
n distinct eigenvalues, there exist a unique directed graph G and a unique 
G -canonical matrix pair (M, A^)can such that (M, A^) is similar to (M, A^)can- 
Thus, (M, A^) is similar to {M', N') if and only if (M, iV)can = (M', A^')can- 
(b) Each pair {M, N) of n x n matrices in which M has n distinct eigen- 
values is similar to a direct sum of indecomposable canonical matrix pairs. 
This direct sum is uniquely determined by (M, A^), up to permutation of sum- 
mands. 



The statement (a) of Theorem 16. II follows from fllSp and f ll9p . The state- 
ment (b) is a consequence of the following remark. 
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Remark 6.1. The direct sum in Theorem 16.1( b) is permutationally similar 
to (M, A^)can and can be obtained from it as follows: The directed graph G is a 
disjoint union of trees; denote them by Gi, . . . , Gg- Let un < Ui2 < ■ ■ ■ < uu^ 
be the vertices of Gi. Let {Ai,Bi) be the pair of ti x ti submatrices of the 
matrices in (M, A^)can formed by rows un, . . . , uu- and columns un, . . . , Uit^. 
Definition 16 . 1 1 ensures that the Uii,Ujk entries of the matrices in (M, iV)can are 
zero if z 7^ j. Therefore, (M, A^)can is permutationally similar to 

(Ai, 500(^2, 52)0- ••e(A„ 5,), 

which is the desired direct sum. Each {Ai, Bi) is a G --canonical matrix pair, 
in which G- is the tree obtained from Gi by relabeling its vertices un, . . . , un^ 
as 1 , . . . , t j . 
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